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We present measurements of the angular momentum flux (torque) in Taylor-Couette flow of 
water between independently rotating cylinders for all regions of the (fli, fta) parameter space at 
high Reynolds numbers, where (0,% ) is the inner (outer) cylinder angular velocity. We find that 
the Rossby number Ro = (fii — Q2) /O2 fully determines the state and torque G as compared to 
G(Ro = 00) = Goo- The ratio G/Goc is a linear function of Ro -1 in four sections of the parameter 
space. For flows with radially-increasing angular momentum, our measured torques greatly exceed 
those of previous experiments [Ji et al, Nature, 444, 343 (2006)], but agree with the analysis of 
Richard and Zahn [Astron. Astrophys., 347, 734 (1999)]. 



Rapidly-rotating shear flows are ubiquitous in geophys- 
ical and astrophysical settings such as planetary atmo- 
spheres, stellar interiors and accretion disks. In order 
for essential fundamental processes to occur, like matter 
inflow towards compact objects pQ, there must be an ex- 
change of angular momentum through such flows. Deter- 
mining the flux of angular momentum in rotating shear 
flows is difficult and has been actively studied pifTTJ) since 
the initial measurements of Wendt [IT] and Taylor [12] . 
While shear tends to destabilize fluid flows, rotation sta- 
bilizes them whenever angular momentum increases ra- 
dially outward (dL/dr > 0, the Rayleigh criterion [13]); 
although, shear turbulence has been shown to occur even 
in Rayleigh-stable systems [14] . Most astrophysical flows 
are Rayleigh-stable, which in the absence of other insta- 
bilities cannot transport angular momentum efficiently 
enough to allow for the breadth of observed phenomena. 

The question remains as to how astrophysical objects 
transport angular momentum at the observed rates [T] . It 
is widely accepted that turbulence is mediated by an in- 
stability, which enhances angular momentum transport; 
although, the particular driving mechanism remains con- 
troversial. In the case of electrically-conducting accretion 
disks, magnetic fields can give rise to the linear magneto- 
rotational instability (MRI) [TShTS] . The MRI increases 
the flux of angular momentum, but it is still unknown 
if (i) it is the only instability and (ii) it is sufficient to 
produce the observed angular momentum transport. 

Ji et al. [5] measured the flux of angular momentum 
in Taylor-Couette flow between two cylinders that ro- 
tated with angular velocities ill and Q 2 for the inner and 
outer boundaries, respectively. They focussed on "quasi- 
Keplerian" flows in water, where the angular velocity 
decreases radially (dfl/dr < 0) while the angular mo- 
mentum increases (dL/dr > 0), as in many astrophysical 
flows. Their measurements of the flux of angular momen- 
tum, inferred from velocity measurements, were indistin- 
guishable from those in solid-body rotation (f2i = £l 2 ), 
where the flux is zero. Ji et al. concluded that there 
was no hydrodynamic instability and that non-magnetic, 



quasi-Keplerian flows at Reynolds numbers up to 2 x 10 6 
are "essentially steady." These results are at odds with 
the analysis of Dubrelle et al. [7] , who used prior velocity 
measurements to infer that the flux of angular momen- 
tum is nonzero for Rayleigh-stable flows. However, there 
has yet to be a study that directly measures the flux of 
angular momentum in all portions of the (£1%, O2) param- 
eter space. Furthermore, Ji et al. could not quantify the 
flux for quasi-Keplerian flows since the random errors in 
their measurements exceeded the measured values. 

In this Letter, we characterize the flux of angular mo- 
mentum (torque) in Taylor-Couette flow between inde- 
pendently rotating cylinders for all regions of the param- 
eter space up to Reynolds numbers Re = (fli — il 2 ) (b — 
a)(a + b)/2v = 4.4 x 10 6 , where v is the kinematic vis- 
cosity and a (b) is the inner (outer) radius. Rather than 
inferring the torque from velocity measurements, as in 
[8], we directly measure the torque required to drive the 
inner cylinder [2]. The precision of our measurements 
allows us to distinguish between solid-body rotation and 
other flow states observed, in contrast to ref. [8]. 

We partition the experimental parameter space (Fig. [I]) 
into four regions. The regions are defined by the Rossby 
number, which we claim is the main controlling parame- 
ter for the dynamics: 

Ro={n 1 -n 2 )/il 2 = {Re 1 /r ! Re 2 )-l, (1) 

where r\ = a/b, Re\ — Q,\a(b — a)/v and Re 2 = fl 2 b(b — 
a)/v are the radius ratio and inner and outer cylinder 
Reynolds numbers. The boundaries co-rotate in region 
I with < f^i < ft 2 . Region II also has co-rotation 
with radially increasing (decreasing) angular momentum 
L (angular velocity J7), as in quasi-Keplerian flows. Re- 
gions II and III are divided by the Rayleigh stability crite- 
rion dL/dr — 0. Region IV is characterized by counter- 
rotation with —4 < Ro < —1. Regions I and II are 
Rayleigh-stable and are separated at Ro = (solid-body 
rotation). Regions III and IV are linearly unstable since 
dL/dr < 0. The division between regions III and IV is 
a result of the observation of a maximum in the mea- 
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FIG. 1. Our experiments span a large range of the (Re2, Re\) 
parameter space, which we divide into regions. Region I (blue 
diamonds) is defined as — 1 < Ro < 0. Region II (red squares) 
has < -Ro < 77 _ 2 — 1 , where Ro — i]~ 2 — l defines the Rayleigh 
stability criterion [T3]- Region III (black circles) is for Ro < 
—4 and n~ 2 — 1 < Ro. Finally, region IV (green triangles) has 
—4 < Ro < — 1. Data are not acquired very near Ro = 0, since 
the torques are comparable to our measurement precision of 
0.01 Nm (G ~ 10 8 at 50 °C). 



sured torque at Ro = —4. The dynamics are symmetric 
under negating both Cli and therefore, regions I-IV 
describe the entire Taylor-Couette parameter space. 

Our experiments use water as the working fluid and are 
conducted in the apparatus constructed by Lathrop et al. 
0, which has been modified to allow the outer cylinder 
to rotate independently. The acrylic outer cylinder used 
in previous experiments [2J |3] has been replaced by an 
anodizcd aluminum cylinder with the same inner radius 
b = 22.085 cm and length L = 69.50 cm. The inner 
cylinder is stainless steel with a radius a = 16.000 cm 
yielding a radius ratio rj = a/b = 0.7245 and an aspect 
ratio r = L/(b — a) = 11.47. The inner cylinder is ro- 
tated up to Qi/2tt = 20 s _1 while the outer cylinder may 
be rotated in either direction up to |fi 2 /27r| = 10 s _1 . 
Both angular velocities are measured precisely by shaft 
encoders and controlled to within 0.2% of the set value. 

The axial boundaries rotate with the outer cylinder, 
although the ideal Couette geometry is unbounded axi- 
ally. To avoid end effects in our torque measurements, 
the inner cylinder is divided axially into three sections of 
length 15.69, 40.64 and 15.69 cm (see Fig. 3 of 2 ). Only 
the central section of the inner cylinder senses the torque 
of the fluid as described in ref. [2]. Therefore 2.58(6— a) 
from each of the axial boundaries, where secondary cir- 
culation setup by finite boundaries (Ekman pumping) is 
strongest, are avoided in the torque measurements. 



The local wall shear stress is measured at the outer 
boundary using a flush-mounted hot film probe. The 
probe is located at the mid-height of the experiment. 
The measurements are performed in the constant tem- 
perature mode using a Dantec mini-CTA anemometer. 
The frequency response of the probe exceeds our sam- 
pling rate of 10 kHz. The shear stress measurements are 
calibrated in situ using the method described in ref. [2]. 

The desired accuracy of our measurements requires 
that the temperature of the water be precisely con- 
trolled. In contrast to prior experiments [3] where 
the system was cooled at the axial boundaries, we con- 
trol the temperature through the outer cylinder. This 
procedure is superior owing to the 6.5 fold increase in 
temperature-controlled surface area. Furthermore, the 
working fluid is now temperature-controlled along the 
entire axial length of the experiment. This is particu- 
larly important for the flows in regions I and II of the 
parameter space, where mixing is greatly reduced. In 
all of our measurements the temperature is controlled to 
within 0.02 °C of 50 °C yielding a kinematic fluid vis- 
cosity of v = 0.0055 cm 2 /s, except for Re > 2 x 10 6 
where T = 90 °C and v = 0.0032 cm 2 /s. This control 
algorithm and temperature range would not be possible 
with an acrylic outer cylinder owing to the poor thermal 
properties compared to those of aluminum. 

We study the scaling of the torque as a function of Re 
for several values of Ro. The measured torque r is made 
dimensionless by defining G — rj pis 2 L c , where p is the 
fluid density and L c the length of the torque-sensing cen- 
tral section of the inner cylinder. Our measured values 
of G are shown in Fig. [2]ja). We note that both Re and 
G are negative in region I and we therefore plot their 
absolute values. We compare our data to the best fit of 
G(Ro — oo) (solid curve) given in ref. [3], which is well 
described by the following: 

Re^ , 

-=L = 1.56 log VG»- 1-83. (2) 

V ""oo 

Our measurements of G obey the scaling observed in pre- 
vious experiments for Ro — oo 13] , even above the pre- 
vious maximum of Re = 1.2 x 10 6 . However, depending 
upon Ro, the value of G for a given Re may be higher or 
lower than G^ [71 [TO] . The main observed dependence 
on Ro is a vertical shift in this representation. 

To determine this shift in G with Ro we measure 
G/Goo, as in ref. [Jj, as a function of Ro' 1 , where Goo 
is given in eq. ^ [3]. Figure [2jb) shows that G/Goo is 
essentially constant for each value of Ro. Therefore, the 
value of Ro fully determines the basic state of the flow, 
which then scales with Re in the same manner as the 
case of outer-stationary Taylor-Couette flow {Ro = oo). 

The behavior of G/Goo is distinct in the four regions of 
the parameter space. For Rayleigh-stable flows (regions 
I and II), G/Goo is less than 0.22. On the other hand, 
for —0.37 < Ro' 1 < the torque is enhanced. Since 
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FIG. 2. Experimental measurements of the dimensionless 
torque G as a function of (a) Re and (b) 7?o _1 with the sym- 
bols defined in Fig. [7J The solid line in (a) represents the fit 
Goo for Ro — co given in eq. |2| [3[. The solid lines in (b) 
correspond to the fits given in eqs. |3]h[6|. 



G/Goo scales as Ro^ 1 within each region, we perform 
linear regressions of G/Goo as a function of Ro^ 1 . The 
resulting functions are shown as solid lines in Fig. [2jb) : 



region I 
region II 
region III 
region IV 



G/Goo = 0.078RO- 1 +0.27 
G/Goo =0.17 ±0.05 
G/G^ = -0.95.Rcr 1 + 1.00 
G/Goo = 2.24RO- 1 + 1.83 



(3) 
(4) 
(5) 
(6) 



co/Q 




In addition to the distinct torque scaling in regions I— 



FIG. 3. Fluctuations in the wall shear-stress strongly depend 
upon Ro. The logarithmic spectral power density is indi- 
cated by color and the spectral frequency ui is normalized by 
Q,2/2n = 8.22 s _1 . The experimentally inaccessible range of 
— 1.25 < Ro < 0.75 is shown as white. 



IV, the wall shear-stress spectra also show marked 
changes with Ro. Figure [3] shows a spectrogram for —2 < 
Ro < 2.1. For Rayleigh-stable flows (-1 < Ro < 0.905), 
the system is characterized by narrow-band, weak shear 
stress fluctuations. The fluctuations are stronger in re- 
gion II, which also shows strong wave modes near Ro = 
0.5. Region III only has strong, broadband fluctuations 
for Ro > 0.95, even though the system becomes linearly 
unstable at Ro — 0.905 in the case of vanishing viscos- 
ity. Finally, the spectra for flows in region IV are also 
broadband but with much stronger wave modes evident. 

Our measurements of the dimensionless torque G and 
wall shear stress spectra indicate that the dominant con- 
trol parameter for rotating shear flows is the Rossby num- 
ber Ro. We have not observed any transitions or evidence 
for nonlinear instabilities with increasing Re in regions I, 
III or IV, although we have indications of hysteresis in 
region II. This likely indicates that quasi-Keplerian flows 
can be nonlinearly unstable, but more systematic studies 
are needed to determine if this is indeed the case. 

In ref. [9] the variations of the torque with rotation 
rates are modeled using exact relations derived from 
the Navier-Stokes equations and assumptions about the 
torque contributions from the radial and vertical veloc- 
ity fields (the "wind"). Our observed variations with Ro 
suggest that, while the exact relations in ref. [9 contain 
only fii — O2, there is a dependence of the wind on Ro. 
The data presented here may be used to determine this 
dependence and provide testable predictions for the ve- 
locity fields (which are not presently accessible). 

Ji et al. [5] did not directly measure G, but instead 
used their velocity measurements to determine the pa- 
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FIG. 4. Scaling of ft, given in eq. Q, and the skin friction 
coefficient c/ = G/Re 2 with i?e. The symbols correspond to 
the portions of parameter space defined in Fig. [T] 



rameter ft, which has been used to interpret angular 
momentum transport in astrophysical objects [2HZ]- In 
this prescription ft = vg/v r f / (vg 2 q 2 ) , where overbars 
indicate mean quantities, vg is the azimuthal velocity, 
vgi = vg — vg, v r i = v r — v^ and q — — <91nf2/<91nr. This 
model describes angular momentum transport as a diffu- 
sive process with a diffusivity ^turb = ft \r 3 dVl/dr\ . Thus, 
for a given f2(r), ft determines the transport of angu- 
lar momentum, with larger values of ft corresponding to 
larger fluxes of L. For quasi-Keplerian flows (region II), 
Ji et al. measured ft = (0.72 ±2.7) x 10~ 6 . This value 
is smaller than the value of ft = (1.5 ± 0.5) x 10~ 5 deter- 
mined by Richard and Zahn [4] for Ro = — 1 (£!i =0). Ji 
et al. attributed the disparity to Ekman circulation pro- 
duced at the axial boundaries, which is reduced in their 
experiments by independently controlling the angular ve- 
locity of the axial boundary 19]. 

To determine ft from our measurements of G we use 
the expression given by Dubrulle et al. [7J 
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where Rc = 2(6 — a)/{b + a) and we take S\ &m /S = 3 
from [7J. Our measured values of ft are shown in Fig. 
[4j We note that ft is proportional to the skin friction 
coefficient Cf = G/Re 2 , which is denoted on the right 
vertical axis. The values of ft are shifted vertically for 
a given Ro by the same factor as G/G^ shown in Fig. 
[2](b). For Ro = -1 we measure ft = (1.84 ±0.03) x 10~ 5 , 
which agrees with the value determined in ref. For 



the flows in region II we determine an average value of 
P = (1.7 ± 0.2) x 10~ 5 , which is markedly higher than 
the value given in ref. [8]. For Rayleigh- unstable flows, 
though, Ji et al. report j3 > 10 -3 [S] whereas our values 
span 2 x 10~ 5 < j3 < 2 x 10~ 4 . 

In conclusion, we have presented the first characteriza- 
tion of the flux of angular momentum (torque) between 
independently rotating cylinders for all regions of param- 
eter space. The reduction or enhancement of the torque 
G/Gqo at a given Reynolds number only depends upon 
the Rossby number Ro. The i?o-dependence of G/Goo is 
well described by eqns. ([3])-(|6|. In contrast to ref. jH] but 
in agreement with ref. [4]7 our measurements of ft, which 
may be used to model angular momentum transport, are 
nonzero for Rayleigh-stable flows. This disparity likely 
indicates that multiple states are possible for Rayleigh- 
stable flows, with our measurements representing a "tur- 
bulent state" and those in ref. [8] a "laminar state." This 
is particularly important for astrophysical flows where 
such nonlinear instability could explain the observed an- 
gular momentum transport. Systematically perturbing 
Rayleigh-stable flows while measuring the torque could 
be used to directly test for nonlinear instabilities. 
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